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Abstract 

In this work, we investigate gradient coupling effect on amplitude death in an array of N cou- 
pled nonidentical oscillators with no-flux boundary conditions and periodic boundary conditions 
respectively. We find that the effects of gradient coupling on amplitude death in diffusive coupled 
nonidentical oscillators is quite different between those two boundaries conditions. With no-flux 
boundary conditions, there is a system size related critical gradient coupling within which the 
gradient coupling tends to monotonically enlarge the amplitude death domain in the parameter 
space. With the periodical boundary conditions, there is an optimal gradient coupling constant ro 
to realize largest AD domain. The gradient coupling first enlarges then decreases the amplitude 
death domain of diffusive coupled oscillators. The amplitude death domain of parameter space are 
analytically predicted for small number of gradient coupled oscillators. 

PACS numbers: 05.45.Xt, 05.45.-a 
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I. Introduction 

The model of coupled nonlinear oscillators provides a simple but powerful paradigm for 
understanding of collective behaviors such as emergent behavior which are widely explored in 
the interacting large number of natural oscillators. Therefore, the study of coupled nonlinear 
oscillators has become a very hot topic in nonlinear sciences and many other interdisciplinary 
fields such as physical, chemical, biological, and even social sciences[l-3|. 

Ensembles consisted of different types of coupled oscillators exhibit various of collective 
behaviors. Among, the amplitude death(AD), which refers to a situation where individual 
oscillators cease oscillating when coupled and go to an equilibrium solution instead, has 
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been actively investigated since the appearance of Ref j4| . AD plays a crucial role in a 
real systems, for instance, it has been extensively found in synthetic genetic networks 
where the AD implies a constant protein expression and the system multi-stability with AD 
is believed to improve the adaptability and robustness of the cellular population. Significant 
progress was achieved in theoretical and numerical analysis of AD in systems of oscillators 
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with various coupling schemes such as all-to-all coupling [8], diffusively coupling [qHuJ even 



in the complex networks where the effects of topological properties of the network on partial 
AD dynamics are explored [12I. AD can be eliminated by introducing random links 

and the influence of spatial disorder on AD in oscillator arrays with local couplings was 
extended. The desynchronizat ion- induced AD are weakened considerably by introducing 
the random deviation from a linear trend of frequencies in array of diffusively coupled limit- 
cycle oscillators with a regular monotonic trend of natural frequenciesjol. 

The gradient coupling, one of anisotropic coupling, is of practical importance in many 
situations, such as in hydrodynamics flows with sloping channels and in plasma systems with 
electromagnetic fields. The gradient coupling has effects in the control of spatiotemporal 
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chaotic systems Q,[l5| and in the synchronization of coupled nonlinear oscillators 
is observed that as the gradient coupling constant increases, the networks' synchronizability 
is possibly enhanced or decreased or even a optimal value of coupling constant exists for best 
synchronizability 19l-l2l|. Therefore, it is also significant to explore the effects of gradient 
coupling on another collective behavior, the AD dynamics, in array of diffusively coupled 
oscillators. However, most of previous studies on AD have been confined only to the cases of 
coupled oscillators with a homogeneous identical coupling from their neighbors, and the effect 
of anisotropic coupling has seldom been studied until in the recent published article |22|, the 
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effects of gradient coupling on the time-delay induced AD is explored. The gradient coupling 
tends to monotonically reduce the domain of delay-induced death island. Typically for the 
occurrence of AD, one of the following two conditions is needed: the time-delayed coupling 



23|or the parameter mismatches lOl. Il2l . |l3|. In addition to these two general conditions, 



recent studies revealed that AD may also happen by dynamic coupling [2J] or conjugate 



coupling [251, |26|]. Since frequency mismatches is widely existing in the natural world. It 
is meaningful to ask how does the gradient coupling influence the frequency mismatches 
induced AD in the array of diffusively couple oscillators. Is the anisotropic coupling beneficial 
to decrease the minimal frequency mismatch needed for AD for given diffusive coupling? To 
answer those question, the influences of gradient coupling on the AD domain of oscillators 
with frequency mismatches are explored. We find that the effects of gradient coupling on 
AD in diffusive coupled oscillators is strongly related to the boundary conditions which has 
effects on the synchronization ability of the diffusively coupled oscillators as discussed in 
Ref . j27l| . With no- flux boundary conditions, the gradient coupling constant monotonically 
enlarge the AD domain until it is larger than the system size related critical value Tc. with 
periodic boundary conditions, the increasing gradient coupling first enlarge then decrease 
the domain of AD. There is an optimal system size related gradient coupling constant to 
realize largest AD domain. The parameter space of AD are analytically predicted for small 
number of gradient coupled oscillators. 
II. Gradient coupling model 

Consider the general form of N gradient coupled nonidentical oscillators with 

X, = /(a;^,X,) + (e + r)(X,+i-X,) + (e-r)(X,_i-X,) (1) 

for j = 1, ...,N{N > 3), where Xj represents the state vector of the jth element; e and r 
are the diffusive and gradient coupling strengths respectively. The uncoupled units Xj = 
f{u)j,Xj) have non- stationary behaviors and meanwhile, accompanying an unstable focus 
X*. Uj is the natural frequencies of uncoupled oscillators. With the frequency mismatches, 
the AD may occur in certain range of coupling constant region when the coupling interaction 
turns the formerly unstable focus X* stable. The distribution of the frequencies has great 
effect on the AD dynamics, i.e. Ref. qI, the AD in the diffusively coupled oscillators 
with a regular monotonic trend of natural frequencies can be weakened by the disorder 
of frequencies distribution. The transition process from partial AD to complete AD in 



3 



diffusively coupled oscillators with frequency mismatches are explored in Ref. 1^. For 
the random distribution frequencies, the frequency deviation has greatly influenced the AD 
behaviors. Without losing generality, we consider coupled nonidentical Landau-Stuart 
oscillators with the no- flux boundary condition z^^i = z^iZq = Zi and the periodical 
boundary condition zjq+i = zi,zq = zjq respectively. The coupled system is presented as 
follows. 



Zj = (1 + ii^j + \zj\ )zj + (e + r){zj+i - Zj) + (e - r)(zj_i - Zj),j = 1, 



(2) 



where i is the imaginary, zj is the complex variables, Wj are the intrinsic frequencies of 
single uncoupled oscillators. For simplicity, we suppose the coupled oscillators have a regular 
monotonic trend of natural frequencies Wj = Wi + (j — j = 1, 2, A^, wi = 1,6ijJ > 0. 
In the absence of coupling, each oscillator has an unstable focus at the origin \zj\ = and 
an attracting limit cycle Zj{t) = e*"^^* = x{t) +iy{t) with different oscillating frequencies Uj. 

For r = 0, Eqj2] is diffusive coupled system as discussed in Ref. jol, Q, where the 
necessary condition of AD is e > 0.5. To explore the influence of gradient coupling on the 
domain of AD, the stability of the complete AD can be analyzed by linearizing the Eq. |2]at 



0,J 



1,2,..., AT. 



rij{t) = (1 - 2e + tUj)r],it) + (e + r)r]j+,{t) + (e - r)r],_,it) 



(3) 



With denotations of ri{t) = {rii{t), r/2(t), Vnit))' , the above equations can be rewritten as 
follows. 



(4) 



where for non-flux boundary condition 
/ 1 — 2e + ioji e + r 



B = Br, 



\ 



\ 



e — r 1 — 2e + iuj2 e + r, 

e — r 1 — 2e + iu^ 

e — r 



e + r 



1 —2e + iujn j 



or for periodical boundary condition 
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B = Br, 



^ 1 — 2e + iui e + r e — r 

e — r 1 — 2e + iu2 e + r 

e — r 1 — 2e + ico^ e + r 



\ e + r e — r 1 — 2e + icuN / 

Assume that B can be diagonalized by the matrix P: 

p-'BP = dtagi\o,\i,-AN-i) (5) 

where A^, = 0, 1, — 1 are the eigenvalues of the matrix B. The necessary condition for 
AD of Eq|2]is that all real part of the eigenvalues Re{\k) <0, /c = 0, 1,...,A^— 1. Therefore, 
the region of AD state is completely determined by the critical line of all Re{Xk) < 0,k = 
0, 1, A^ — 1. However, it is difficult to diagonalize the matrix B analytically when A^ is 
large. To investigate the effects of gradient coupling on AD domain, we mainly resort to 
numerical simulations for large A^. Moreover, the AD domain can be analytically presented 
for small system size A^. 

III. The non-flux boundary condition 

Let's firstly consider the gradient coupling effects on AD in the coupled system with 
non-flux boundary condition. The AD domain of parameter space e ^ 6uj are calculated 
for different A^ , r in Eql2] as shown in Fig{T](a)~(d) for (A^ = 3,4,10,100) respectively. 
The AD domain is the right and above part of area enclosed by e > 0.5 (red solid line) 
and the curve line of corresponding r. The area enclosed by e > 0.5 and the black solid 
curve for r = (marked with AD) in Fig{T] are the AD domain for r = which is right 
the AD domain of diffusive coupled system. Interestingly, we flnd that the AD region is 
monotonically expanding as the gradient coupling r is gradually increasing until r = r^, 
where Vc is related to the systems size A^. Thus, for an arbitrary given system size A^, 
the gradient coupling may monotonically increase the AD domain for r < Tc. However, 
when r > Tc, the AD domain becomes the area enclosed by all e > 0.5 and 6u) > and 
keeps constant for increasing gradient coupling strength r. Therefore, the gradient coupling 
tends to minimize the frequency mismatch needed for AD in diffusively coupled oscillators. 
According to FigH], there is another remarkable thing that smaller frequency mismatches 600 
and larger diffusive coupling strength is needed for AD in the gradient coupled system with 
larger system size A^ for arbitrary given constant r. However, based on these observations. 
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we may predict the critical gradient coupling rc for each A^. A normalized scaling factor 
is defined as R{r) = 1 — S{r)/S{0), where S{r) denotes the area of non-amplitude death 
island in the parameter space e ~ for e > 0.5. Obviously, -R(O) = 0. The relationship 
between R{r) and r are presented in FigJ2]^a) for different system size A^. R{r) monotonically 
increases with increasing r until r = r^, R{r^ = 1. The critical gradient coupling has 
positive linear relationship with the system size A^ as shown in Figj2]^b). 

The AD domain of the gradient coupled system can be analytically presented according 
to Re{Xi) < 0(the eigenvalues of matrix Bn) when the size A^ is small, for example, A^ = 3, 4. 
When A^ = 3, the eigenvalues of matrix Bn are presented as follows, 

Ai = 1 -2e + i{l + Suj), 
A2,3 = 1 - 2e ± V2e2 - 2r^ - 6oj^ + t{l + 600) (6) 

Let Re{Xi) < 0, the boundaries of the AD domain can be determined by following equations. 

(1) . Area I, if 

5u'' > 2^ - 2r2 (7) 

then e > 0.5 

(2) . Area II, if 

5u'^ < 2e^ - 2r2 (8) 

then 

e > 0.5, 

600^ > I- 2(e - If - 2r^ (9) 

Therefore, the AD domain for arbitrary given r is composed of area I and II(AD domain for 
r = 0.3 is presented in Figj3]^a), which is completely consistent with the numerical results 
as presented in Figj3t^b)). As r increase from zero to the AD area II enclosed by Eq|9] 
enlarges monotonically. As r > rc, AD is stable for all e > 0.5 and 5uj > Q and keeps constant 
to increasing r. From EqlHl the critical gradient coupling rc of A^ = 3 can be theoretically 
calculated as du"^ > for all e > 0.5, that is, 1 - 2r2 = 0,(rc = ^/T/2 = 0.707). 
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In the case of = 4, the eigenvalues of matrix Bn can be given as 
Ai,2,3,4 = l-2e±0.5v^ + ^(l + 35a;/2), 

The boundaries of the AD domain can also be presented I, II, and III according to 

i?e(Ai,2,3,4) < 0. 

(1) Area I, enclosed by EqJTT] 

e > 0.5 (11) 

(2) Area II, enclosed by EqiH 

e > 0.5, 

Stu^ > -(-20 - + 80e - 746^ 
9 

+2Vl584e2 - 512e - 2240e3 + 192er^ + 64 - 48r^ + 1189e4 - ISSeV^ - 27r^) (12) 

(3) Area III, enclosed by EqJT3l 

e > 0.5, 

(e^ - r^)/2 < 6uj^ < (e^ - r2)/0.4, 

6u^ > -(-10-6r2-34e2 + 40e 
4 

+2V4r4 + 16e2r2 - 24er2 + 6r^ + 1246^ - 264e3 + 210e2 - 72e + 9 (13) 

The AD domain for r = 0.5 is marked I, II and III in FigOJ^c) which is coincided 

with the numerical results in Figj3]^d) well. According to EqJT2l the critical gradient coupling 
Tc of = 4 can be determined by EqUH 

-(-20-6r2 + 80e-74e2 
9 

+2Vl584e2 - 512e - 2240e3 + 192er2 + 64 - 48r2 + 1189e4 - 138e2r2 - 27r^) = (14) 
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The solution of Eq{Tl] is EqJTSl 

e = 1 ± — ± + 0.3^5 - (1 + OAVby (15) 

Since e is real, then 0.7 + 0.3v^ - (1 + 0.4v^)r2 > 0, that is r, = -^^^^ = 0.8507. 
IV. The periodical boundary condition 

Now let's explore the gradient coupled oscillators with periodical boundary condition. 
By checking the AD domain of parameter space e ~ (5w which is similar to that in Fig. [1] 
for = 3,4, 10, 100 as shown in Fig|l](a)~(d), where the AD region is the above part of 
area enclosed by line e > 0.5 and the curve line of corresponding r. One may find that when 
A^ = 3,4, the gradient coupling tends to decrease the AD domain monotonically as shown 
in FigJll^a)(b), while for A^ = 10, 100, there is a critical gradient coupling constant r^, within 
which the gradient coupling first increase then decrease the AD domain. That is to say, there 
is an optimal gradient coupling constant with which the coupled system has the largest 
AD domain in parameter space e ~ 6u). However, when r > r^, the gradient coupling begin 
to shrink the AD domain of formally diffusively coupled system. To quantified the critical 
gradient coupling constant Vc and the optimal gradient coupling constant for different 
size A^, a factor is defined as R'{r) = 1 — S{0)/S{r), where S{r) is the area of the non-AD 
domain of corresponding r for e G (0.5, 20), 5a; G (0,20). (Obviously, -R'(O) = 0; ^^^^^^ < 
means the increasing r enlarges AD domain; When R'{r) > -R'(O), the gradient coupling have 
smaller AD domain than the formally diffusively coupled system) . Then Tc is determined by 
R'i^c) = and the optimal gradient coupling constant are determined by ^^^^^^ \ r=ro = 0. 
The relationship between R'{r) and r for different system size A^ are presented in Figj5]^a). 
Obviously, there is system size related and Vc- The detail relationship between and Vc 
and system A^ are presented in Figj5]^b). Vc is linearly increase with A^. = for A^ = 3,4 
and increases nonlinearly with increasing A^. 

Accordingly, the AD domain can be predicted for small system size for example A^ = 4. 
The eigenvalue of Bp can be presented as follows, 

Ai,2,3,4 = l-2e±0.5y^ + i{l + 1.5Suj), 
Ti,2 = 8e2 - 8r2 - 55a;2 ± 4^^, 
S = 6u^- 4Su\e^ - r^) + 4{e^ + r^. (16) 

The boundaries of the AD domain are determined by -Re(Ai_2,3,4) < 0. The value of 
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-Re(Ai, 2,3,4) has various forms for corresponding parameters. 
(1) if T < 0, i.e. 



Su^ > -i-2r^ + 2e^ + 2Vr'^ + 3Ae^r^ + e^) (17) 
9 



then -Re(Ai_2,3,4) = 1 — 2e 
(2) if T > 0, i.e. 



Suj^ < -{-2r^ + 2e^ + 2Vr^ + 346^2 + e^) (18) 
9 



then 



Ti,2 = - Sr^ - 56uj^ ± 4^6uj^ - A6u\e^ - r^) + A{e^ + r^y (19) 

Thus, the critical hne of AD domain can be found by solving the equations -Re(Ai^2.3,4) < 0. 
The AD domain consist of two areas 

(1) Area I enclosed by solution of e > 0.5 and EqlTTt 

(2) Area II enclosed by solution of EqMland Eqj20i;i?e(Ai,2,3,4) < in EqUH]) 

g{e,r) = |(-18e2 + 20e-2r2-5 



+2 V81e4 - 144e3 + 186^2 + 100e2 + IQr^e - 32e + - Ar^ + 4) (20) 

The AD domain of r = 5, 10 is presented 1 and II in Figj6](a)(b), where area I is 

enclosed by e > 0.5 and EqJTTl while area II is enclosed by Eq{T8]and Eqj20l It is also well 
consistent with the numerical results as shown in Figl6t^c)(d) respectively. 

It is necessary to make some discussion. For convenience of analysis, the frequency 



distribution is set as regular monotonic trend as in Ref.^|. To explore the influence of 
frequency distribution on the AD domain, the noise is added to each frequency. Set wj = 
wi + {6u{j — 1)+^), where ^ is the gauss noise with strength a, namely, (^j) = 0,{C,i{t)^j(t')) = 
a6ij6{t — t'). The AD domain of = 10, cr = 2,r = 5 is presented in Figl71^a)(b) for different 
realization of noise. Where the noise may either enlarge or shrink the AD domain for different 
realization of noise. To see the effects of the noise strength a on AD domain area, the factor 
C((j) = 1 — S{(r)/S{0) versus a is presented in FigJTl^c), where S{(t) is the non-AD domain 
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area for noise constant a for e G (0.5, 1.5), G (0,2). Obviously, with increasing noise 
intensity, the noise tends to enlarge the deviation of the factor C(cr), the effects of noise on 
the AD domain deviation is positively related to the noise intensity a. Secondly, we may 
point out that the gradient coupling effects on AD domain is not exclusively existing in the 
coupled periodical oscillators but also existing in gradient coupled chaotic oscillators such as 
Rossler system dx/dt = wy—z, dy/dt = wx+0.165z, dz/dt = 0.2+z(x— 10) with the coupling 
scheme as described in Eq{T|, where wj = l+{j—l)6ijj, j = 1, 2, A^. The factor R'{r) defined 
above {S{r) is the non-AD domain area for corresponding r for e G (0.5, 1.5), 6uj G (0,2)) 
versus gradient coupling constant r for A^ = 10 coupled Rossler oscillators with periodical 
boundary condition is presented in FigJTl^d) where the increasing r firstly decreases then 
increases R'{r), which indicate that the gradient coupling first enlarge then shrink the AD 
domain area of coupled Rossler system. Moreover, there is also an optimal gradient coupling 
constant = 0.1 and a critical gradient coupling constant = 0.124. 
Conclusion 

The effects of gradient coupling effect on AD domain in an array of coupled nonidentical is 
strongly related to the boundary conditions. With no-flux boundary conditions, the gradient 
coupling tends to monotonically enlarge the AD domain in the parameter space e ~ ^cu 
within the critical gradient coupling constant Vc- When r > the gradient coupling has 
no effects on AD domain any more. With the periodical boundary conditions, there is also 
a critical gradient coupling constant r^. When r > Tc, the gradient coupling tends to shrink 
the AD domain. When r < rc, it first enlarges then decreases the AD domain, that is, there 
is an optimal gradient coupling constant Tq to realize largest AD domain for > 5. The 
AD domain of those gradient coupled system are analytically predicted for small system 
size number. The remarkable thing is that the effects of gradient coupling on frequency 
mismatches caused AD is quite different with that on the time-delay induced AD 2^ where 



gradient coupling monotonically decreases the AD domain till completely eliminates AD. 
Ref. 91] point out that the frequency mismatch 5uj is beneficial to the AD, that is, larger 5u: 
is easier to become AD in diffusively coupled array of oscillators. In this work, one may find 
that the gradient coupling is helpful to AD dynamics since smaller frequency mismatches 5u} 
is needed for AD with proper gradient coupling constant. The gradient effects on the AD 
in coupled non-identical oscillators would be practically valuable to the dynamics control. 
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FIG. 1: (Color online) The critical line of the AD in coupled system with non-flux boundary 
condition. The AD domain is the right and up part of the area enclosed by the line e > 0.5, 6uj > 
and the critical lines of corresponding r. (a) = 3,r = 0(black line),r=0.3 (green line),r=0.5 
(blue line),r=0.7 (cyan line); (b) A^ = 4,r = (black line),r=0.5 (green line),r=0.7(blue line),r=0.8 
(cyan line); (c)A^ = 10, r = (black line),r=l (green line),r=1.5 (blue line),r=1.7 (cyan line); 
(d)A^ = 100, r = (black line),r=l (green line),r=2 (blue line),r=3 (cyan line),r=5 (magenta line). 
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FIG. 2: (Color online) (a) The normalized scaling factor R{r) vs r for iV = 3, 10, 50, 80, 100, 200, 
respectively. R{r) = 1 — S{r)/S{0) with S{r) standing for the area of non-AD island in parameter 
space e ^ 6uj for e > 0.5. For each given system size N, a monotonic increasing of R(r) vs r 
and a critical value Tc are clear. All the data are numerically obtained by directly integrating the 
coupled system Eql2j (b) The critical gradient coupling constant rc vs system size N. Tc is linearly 
increasing with N. 
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FIG. 3: (Color online) The boundary line of the AD domain in coupled system with non-flux 
boundary condition, (a) Analytical results of the AD domain for N = 3r = 0.3(area I and II). 
magenta dashed line is EqjTl green dashed line is EqjUJ (b) Simulation results with parameters 
corresponding to (a), (c) Analytical results of the AD domain for = 4, r = 0.5(area I, II and III), 
(d) Simulation results with parameters corresponding to (c) 
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FIG. 4: (Color online) The critical line of the AD in coupled system with periodical boundary 
condition. The AD domain is the right and up part of the area enclosed by the line e > 0.5 and 
the critical curve lines of corresponding r.(a) r = 0, 5, 10, 15 ,N = 3. (b) r = 0, 5, 10, 15, N=4. (c) 
r = 0, 0.4, 1, 5, 10, N=10; (d) r = 0, 0.5, 1, 2, 5, 10, 15,N=100. 
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FIG. 5: (Color online)(a) The factor R'{r) vs r for iV = 10, 20, 50, 80, 100, 200 respectively. R'{r) = 
1 — S{0)/S{r) with S{r) standing for the area of non-AD island in parameter space e ^ 6uj for e G 
(0.5, 20), (^w G (0, 20). For each given system size N{N > 5), R'{r) firstly decrease to negative then 
increase to positive. There is a optimal gradient coupling constant To which satisfy R'{ro) = R'min 
and a critical gradient coupling constant Tc which satisfy R'{rc) = 0. (b) The optimal gradient 
coupling constant Tq and the critical gradient coupling constant Vc versus system size N{ro = 
when iV = 3,4). 
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FIG. 6: (Color online) (a) (b) Analytical result of the critical line of the AD domain in coupled 
system with periodical boundary condition and = 4. r = 5,r = 10 respectively. The AD domain 
is consisted of area I and area II. (c)(d) The simulation results of AD domain with parameter 
corresponding to (a) (b) respectively. The analytical results are consistent with the simulation 
results. 
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FIG. 7: (Color online) (a)(b) The effects of noise intensity on AD domain for gradient coupled 
oscillators with periodical boundary conditions(A^ = 10, r = 5), AD domain of two realizations 
of noise with intensity a = 2 (cyan region). The red curve is the boundary of AD without noise. 
The noise can either enlarge or shrink AD domain. (c)Factors C{a) = 1 — S{a)/S{0) versus noise 
intensity a. S{(t) is the area of non-AD domain with noise intensity a. The deviation of C{a) 
is increasing with deviation of AD domain area. The stronger noise increase the deviation of AD 
domain, (d) R'{r) versus r of the gradient coupled Rossler oscillators with periodical boundary 
conditions and N = 10, where rc = 0.124, = 0.1. 
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